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FIRST-ORDER THEORY FOR UNSTEADY MOTION OF 
THIN WINGS AT SUPERSONIC SPEEDS 
By Barry Moskcwitz and W. E. Moeckel 


SUMMARY 

An approximate method is presented for determining the pres- 
sure distribution due to unsteady motion for thin wings of fairly 
general plan form. The local source strength, which is propor- 
tional to the normal component of the perturbation velocity, is 
represented by a power series in the wing coordinates, in which 
the coefficients are functions of time . The method is valid for 
arbitrary motions when the second derivative of the perturbation 
velocities with respect to time is not large. 

As examples of the method, the load distributions due to 
oscillations normal to the plane of the wing and pitching oscil- 
lations about a spanwise axis are evaluated as functions of time 
for a swept wing with straight supersonic leading and trailing 
edges. For the frequency of oscillation chosen in these examples, 
the magnitude of the loading differed only slightly from that 
obtained by neglecting the tine delays, but the positions of the 
constant -pressure lines were noticeably altered. 


INTRODUCTION 

The problem of predicting the pressure distribution over thin 
wings undergoing unsteady motion at supersonic speeds has been 
studied by a number of investigators. In reference 1, equations 
are given for the velocity potential in the vicinity of thin wings 
for which the top and bottom surfaces are mutually independent. 

An extension to finite wings with interacting top and bottom sur- 
faces is given in reference 2, but the area-cancellation technique 
required therein is shown to be valid only when the second deriva- 
tive of the perturbation velocities with respect to time is not 
too great. 

For such relatively slow motions, however, a simpler method 
of analysis than that presented in reference 2 can be developed. 
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This method, which is presented herein, consists in expressing the 
local source strength in a power series of the coordinates, so that 
integrals similar to those encountered for steady motion are 
obtained. This method is shown to he accurate to the degree 
of approximation as that given in reference 2 . The a nal ysis was 
completed at the NACA Lewis laboratory during June 1949. 


SYMBOIS 

The following symbols sire used in this report* 
A amplitude of oscillation 

a,a',d,d',d",e integration limits 
functions of time 
wing choid 

pressure coefficient 
velocity of sound 
functions of x, y, and t 


a ,b 
n* n 


k 


slope of straight leading edge in oblique 
coordinates 

Mach number 

area included in forward Mach cone from (x,y) 
time 

free -stream velocity 

oblique coordinates parallel to Mach lines 
oblique coordinates of wing vertex 
oblique coordinates of point on wing surface 
component of perturbation velocity in z -direction 


Cartesian coordinates 
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"w 

P 

no 

cr 

a 


a,b 


vertical coordinate of any point on wing 

angle of chordline with respect to free -stream 
direction 


cotangent of Mach angle, \jvP-i 
half -wing span 

effective local source strength, w/u 

source strength at time t-T plus source strength 
at time t-T^ a 


time delay 


, (*-»» , V(*-n 2 -P 2 (y-.i) 

i O T O 


fi 2 c 


P 2 c 


time delay 


, izZM V(x-» 2 -P 8 (T-n)^ 

' n rt 


$ 2 c 


P 2 = 


<p 

(0 


perturhation-velocity potential 
angular velocity of osoillation 


AKA1YSIS 

For thin wings, the equation for the time -dependent velocity 
potential in the plane of the wing, as given in references 1 and 2, 
is 



y a.b 


*(i,y,t) « - i 

2 *^ M(i-|) 2 -P 2 (y-Ti ) 2 


( 1 ) 


where the source -strength distribution is given by 


0 a,b = + 0(1,7!^-^) 


( 2 ) 
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The integration is performed in the z=0 plane over the area 
included in the forward Mach cone from any point (x,y) and may 
include areas off the wing surface. The time delays T and T, 
are given hy (reference 1) ° 



The functional dependence of a a on these time delays can 
he rendered explicit hy expanding the two terms of equation (2) 
in Taylor's series about T a and Tjj, respectively: 



T a +T b 

2 bt 


1 I r a 2+T h 2 ^ ^ 0 (?,T],t) 

27 \ 2 J Zt 2 

where cTq( | , "n,t ) is the source strength corresponding to zero 

time delay. This representation of the actual source -strength 
distribution O a ^ requires that none of the derivatives of 

aodj^t) are infinite and that the series is convergent. An 

estimate of the magnitude of the angular velocity, for which 
higher-order terms can he neglected, can he obtained hy con- 
sidering an oscillatory motion of the type 




a 0 (|,n,t)co e m 

For such motions, successive derivatives of increase in 

magnitude hy the factor (0. If T a and T ^ are of the order 

M/(1000 p 2 ) seconds (see equation (3)), then successive terms 
of equation (4) decrease approximately in the ratio COM/(1000 P 2 ). 
If an error of approximately 5 percent is considered allowable in 
the representation of O a •j ) , then the neglect of all except the 

zero- and first-order terms may he considered valid for angular 
velocities less than approximately (200 P 2 )/m radians per second. 
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In reference 3, the zero-order approximation for unsteady 
motion is mentioned. The solution for this approximation, as 
seen from equations (4) and (1), is at each Instant equivalent 
to the steady -state solution corresponding to the instantaneous 
distribution of source strengths. The solution for the first- 
order approximation is considered herein. (Special cases for 
this approximation are discussed in references 4 and 5 but the 
derivations are not general.) The source -strength distribution 
for this approximation becomes 


'a,b 


_ , (x-l)M 

- “2 St 


P 2 C 


(5) 


In order to obtain solutions for finite wings with interact- 
ing upper and lower surfaces, the procedures developed by Evvard, 
whereby the effect of disturbed flow outboard of the wing is 
replaced by an equivalent integration over the wing surface, can 
be employed. In reference 2, this area-cancellation technique, 
which was originally developed for steady motion, is Bhcwn, in 
effect, to be valid also for unsteady motion if the motion is 
such that the square -root terms in the time-delay expression 
(equation (3)) are negligible in the equation for C & ■ t) . Inasmuch 

as these square -root terms do not appear in the zero-order and the 
first-order approximations given in equation (5), the area- 
cancellation technique is entirely valid up to the first-order 
approximation. No general procedure is yet available for higher- 
order approximations when the upper and lower surfaces of the wing 
interact . 

Equation (5) leads to an essential simplification of the 
procedure described in reference 2 for types of motion for which 
the area -cancellation technique is valid. This simplification 
consists in expressing the velocity potential in terms of integrals 
commonly encountered for the steady motion of wings. If, for 
example, the motion can be represented as a power series In the 
coordinates, such as 


> * 1 ^) 



M“ + Vi n > 


( 6 ) 


then equation (5) becomes 
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vhere and b n are any function of time consistent with the 

restrictions imposed by retaining only the first two terms of the 
Taylor’s series (equation (4)). Substitution of equation (7) in 
equation (l) results in integrals of the type 



fc r n s dt dri 
x-g) 2 -P 2 (y-il ) 2 


where r and s are integers and the integrals have coefficients 
that are functions of x, y, and t. 


Although motions may occur thaij require values of n greater 
than one in equation (6), many common types of unsteady motion are 
represented by values of n of zero and one. For these motions, 
equation (7) becomes 

Ca,b ■ f (^y^)+s(3C,y,t)|+h(x,y,t)T>fi(i,y,t)5Ti+j(x,y,t)t 2 (8) 


where 


f 


+ b rt - 


jdM 

p2 

P o 


dao 

dt 


S 


a l 


^1 . M ( ^0 

p 2 0 St + p 2 0 "l “5fc + 



h 


2 



P 2 C 
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2M ^1 
p 2 c ^ 

2M ^ a l 

0 2 c & 


If equation (8) is used in equation (1), the quantity f at any 
instant is seen to be an effective angle of attack, and Ug and 
Uh are effective rates of steady pitch and roll about the | =0 
and q=0 axes, respectively. For a ving of the type for vhich 
the methods of reference 2 are applicable, the procedures used to 
determine the velocity potential in steady flight at angle of 
attack and in steady roll and pitch may be used for all unsteady 
motions representable by the first three terms of equation (8) . 
When the velocity potential has been determined, the load dis- 
tribution is obtained from the equation 




1 Sq> 


vhich, from equations (8) and (l) becomes 


C p = “ U 


^2 »4 


cp cp, fk + i sA 

^2 u btj l^dx + U ht/ V4 \9ir XT bbj 


*5 + 


where 





dt dn 




( 9 ) 
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*2 - "ST 


JL / / ^ 

V(x-I) 2 -P 2 ( y-n ) 2 


s 




ti dl dn 

2 * ' \/t=c-t) 2 -P 3 (y-n ) 2 

"S 



_ u 
^4 = - 2^ 



|tj a| dTi 


_ u 

9 5 “ " 2 * 



V(^-E) 2 -e 2 (y-Ti) 2 
I 2 al 

f\/(x-£) 2 -P 2 (y-T}) 2 


The first term within the brackets of equation (9) yields a pressure 
coefficient corresponding to steady flight at angle of attack f; 
the second term is the pressure coefficient due to steady pitch at 
instantaneous rate Ugj and the third term is the pressure coef- 
ficient corresponding to steady roll at instantaneous rate Uh. 
Equations for these pressure coefficients are derived for a fairly 
general class of wings in reference 6. The remaining terms in 
equation (9), however, are not, for the most part, to he found in 
steady -motion literature and must he independently evaluated. 


APPLICATION 

As an illustration of the procedure outlined in the section 
entitled "Analysis," explicit expressions will he derived for the 
load distribution resulting from motions expressible in the form 


a a,b = + e(x>y/t)! ( 10 ) 

A flat-plate wing of the plan foim shown in figure 1 will he con- 
sidered. This wing is similar to that analyzed in reference 6, 
except that no subsonic trailing edges appear. 
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The integrations required to determine load distribution are 
simplified if, as in reference 2, the equations are -transformed 
to the oblique coordinate system whose axes are parallel to the 
Mach lines. The transformation equations are: 


u = $ (l-N 

V 

* n d+Pi) 

\ = | (v+u) 

T) 

\5 (t ' u > 

u w - gfr ( x -Py) 

v w 

= (x+Py) 

1 “ | ^w +u v) 

r 

* (W> 


(11) 


The elementary area is (2P/M 2 )du dv. Xn this system of coordinates, 
for types of motion given by equation (10), equation (9) becomes 


p _ _ _ 

P “ U 


f 


M_ 

2P 


M *i\ 

\<K + +8 2p 




dx U St 


a* u at 


where 


m U 

= “ 2Mrt^ 


** “ 


du dv 


S V (uw-tt) (v w -v) 


(v-t-u) du dv 


2M k /\J(u v -n) (v v -v) 



( 12 ) 


The integration procedure, which employs area-cancellation 
technique for regions influenced by areas off the wing plan form, 
is the same as that described in reference 6, so that only an 
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oatline of the derivation need he repeated herein. As in refer- 
ence 6, the origin of coordinates is placed at one of the junc- 
tures of supersonic and subsonic leading edges. The straight 
supersonic leading edges are defined by u-j_ = -kv and 


U 2 



respectively, and the curved subsonic lead- 


ing edges are defined by u* u 3 (v) and v » v^u), respectively. 

The Mach lines from the vertex and from the junctures of the 
supersonic and subsonic leading edges divide the plan form into 
eight types of region. The integration limits in eaoh region 
are found as described in reference 6. The presentation is 
simplified if the expression for Cp is given for a general 

region of the wing, together with the values of the limits to be 
used for each region. For the general region the quantities 
required in equation (12) become 



“18 


V Cu w -u) (y Jai -«’io< :L -*>) + ~ 


a' 


2 


(13) 
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f- 


—i _ i 


k ] te -1 Huy-n) 

3/ 4k j 'V v^+kn-MtioCl-k) 


a* 


2^+2 v v +T4+a' 




(13a) 


Bq ?L ^ _ _y_ 
2teiv + W 2 *P 




au/'U V W -V4 


av/ Vu w -U3 


k4l 


-j /%-u“ 

tan V Ev™ 


tan 


kCu^-u) 
v w +ku-M^o(l-k) 


(13*) 


a’. 
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vhere the integration limits for each region are given in the follow- 
ing table: 


Wing regions (fig. l) 


Limit 

I 

n 

in 

IT 

T 

VI 

VII 

VIII 

a 

Ujr 

% 

^w 

a 

u 3 

u d 

u 3 

-kv 4 

d 

EH 

-kv v 

u 3 

u 3 

U3 

u 3 

u 3 

ud 

a' 

*w 

Sr 

Uy 



7 4 

k 

>-»•! 

v 4 

k 

Ud 

d' 

U-w 

% 

Sr 

u d 

n 

u d 

U 3 

u 3 

d" 

”w 


Sr 

*w 

Sr 

v 4 

k 

Mq 0 

-^d-k)- 

Zi 

k 

-kv 4 

e 

v w 

v w 

u 3 
- k 

u 3 

"Y 

-ku 3 + 

: Mq 0 (l-k) 

_ U 3 

k 

-ku, + 
Mq 0 (l-k) 

-kU3 + 

Mqo( 1_k ) 


Substitution of equations (13) to (13c), with the appropriate limits, 
in equation (12) gives the pressure coefficient in each region 
corresponding to unsteady motions of the type of equation (10) . 

From equations (13) to (13c) and (12), load distributions 
corresponding to two types of unsteady motion have been computed 
for the wing shown in figure 2. For this wing, which has stream- 
wise tips and supersonic trailing edges parallel to the supersonic 
leading edges, only the first five regions of the wing shown in 
figure 1 are present. 

' i 

The first type of motion considered is an oscillation of the 
wing in the z-direction at constant attitude Oq with respect to 

the free -stream direction. The source strength as a function of 
time is 
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Ob(t) = = A cos 03t + ot^ (14) 

where A is the amplitude of oscillation. The amplitude A is a 
constant herein hut can also he a function of time. Inasmuch as 
the load, distribution due to steady flight at angle of attack is 
known and oan he superimposed on the load distribution due to the 
oscillation, only the first term of equation (14) will he considered, 
that is, the motion is to he defined hy 

0 Q (t) « A cos cot = a Q (t) (15) 


Comparison of equation (15) with equations (6) and (8) shows that 
for this motion, 


ff a t = 2A (cos cot + 2*5“ sin COt ) - sin cot) | (16) 

' \ 0 2 c / \ p 2 c J 

For the top surface of the wing, the required expressions for the 
coefficients and their derivatives heoame 


1 

U 


1 

u 


f = 

-2A f cos cot + 2SJ5? sin tot ) 

V P 2 C J 

(17) 

g = 

A CO sin cot 
p 2 c 

(17a) 

3x 

- m sin cot 

o2_ 

p c 

(l7h) 

II 

2A£0 f , xMCO , A 

-rr- ( sin cot — cos cot } 

Mc \ P 2 C / 

(17c) 

3& 
3x = 

0 

(17d) 

3t 

cos cot 

p2 c 2 

(I7e) 
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Substitution of equations (17) to (l7e) and (13) to (13c) in 
equation (12) yields the desired expression for the load distri- 
bution as function of wing coordinates and time . 

The second type of motion considered is an oscillation of 
the wing about the spanvise axis, | =0. The appropriate expres- 
sion for the source strength can be obtained by considering the 
equation of motion of the z coordinate of the wing. This equation 
is 


z -At cos cot 

V s 


(18) 


At the wing surface the boundary condition is 

1 dz v dz, 


°°~Wz-0 U dt 


w + 1 Sz w 


■5T' r u'St 


(19) 


so that 


a 0 (|,b) = A cos (Ot - sin cot = ag(t) + ai(t)| (20) 

Me 

Inasmuch as the load distribution due to the first term in equa- 
tion (20) is known from the preceding example, only the equation 


a 0 = a-^t) £ 


(20a) 


need be considered here. The total load distribution for the 
pitching oscillation can be obtained by superposition of the 
load distribution due to vertioal oscillation and that obtained 
for equation (20a) . Comparison of equation (20a) with equations 
(6) and (8) shows that 


°a,b 



sin cot + 


XMCO 

□ 2 - 
P c 


cos 



\P 2 c 2 


oos cot ) | 


( 21 ) 


For the top surface of the wing, the required coefficients and 
their derivatives are 



16 


MCA TN 2034 


2A(b 


g > zspL i sin cat + £51“ cos tot ] 

Mc V p 2 c ) 


J = oos 


- 2Aa ^ 
' P 2 C 2 


cos cot 


1- 5g _ 2A<# 4 0S _ ^0) Bln ^ 


U dfc m 2 c 2 


P 2 o 


1 dj. = _ 2AC0 5 
U eft 


M0 2 c 3 


sin cot 


( 22 ) 

(22a) 

(22b) 

(22c) 

(224) 


The remaining coefficients and derivatives in equation (9) are 
zero. 


The unsteady load distributions for the two types of motion 
represented by equations (15) and (20), were computed as func- 
tion of time for the wing shown in figure 3. The Mach number, 
angular* velocity, and velocity of sound were assumed to be A/2 - ^ 

60 radians per second, and 1000 feet per second, respectively. 

For these values, terms containing or fcr could be neglected 
for both types of motion with an error of less than 1 percent 
in the values of C & This error includes the effect of retain- 
ing only the first two terms of the Taylor's series (equation (4)). 

The omission of terms of order C0 2 /c 2 eliminates equations (17e), 
(22a), (22b), (22c), and (224), and the second terms in equations 
(l7o) and (22) from the computation. For the wing analyzed, the 

aspect ratio is LllpJ . a M = /\[2 and the sweepback angle is 26.6° 
b 

(k = l/3). Nondimens ional coordinates (u^/b, v w /b) were used in 
the computation. 
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The results of the computations are shown in figures 4 and 5. 
For vertical oscillations, the loading at t=0 is the steady-state 
lift distribution corresponding to an effective angle of attack 
equal to the amplitude A (equation (16) ) . This load distribution 
is also obtained for the second type of motion when t=0 (equations 
(20) and (21)). For COt = it/4 (figs. 4(b) and 5(a)), the load 
distributions obtained by retaining the first-order terms in the 
source strengths are compared with those obtained with a zero-order 
approximation. The zero -order load distributions, which correspond 
to a neglect of the time delays, are seen to differ considerably 
from those obtained by first-order theory, although for the low- 
angular velocity used the magnitudes of the loadings are almost the ' 
same. For <ht = tf/2 (figs. 4(o) and 5(b)), the loading resulting 
from first-order theory is seen to be small. The loading is zero 
by the zero-order approximation for the vertical oscillation. The 
load distribution for cot = jt,5rt/4, and 3n:/2 are the negative of 
the distributions for cot = 0,ic/4, and it/2, respectively. 


SUMMARY OF THEORY AND RESULTS 

A method has been presented for determining the pressure 
distribution on fairly general classes of thin wings undergoing 
unsteady motion for which the second derivative of the perturba- 
tion velocity, with respect to time, is not large. The method 
consisted in expressing the local time -dependent source strengths 
that represent the wing motion in a power series of the coordinates, 
so that integrals similar to those encountered for steady motion 
were obtained. 

For a rather general wing with no subsonic trailing edges, an 
explicit expression for the pressure coefficient due to unsteady 
motion was obtained as a function, of time. As examples, the load 
distributions for a swept flat -plate wing, with straight supersonic 
leading and trailing edges and a streamwise tip, undergoing oscil- 
lations in a vertical direction and pitching oscillations about a 
spanwise axis, were evaluated. For the frequency of oscillation 
chosen, the magnitude of the loading differed only slightly from 
that obtained by neglecting the time delays, but the positions 
of the lines of constant pressure were noticeably altered. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, August 10, 1949. 
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Figure 4. - Load distribution due to oscillations in vertical direction^ 
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(b) ut s tr/4. 


Figure 4» - Continued. Load distribution due to oscillations in vertical direction- 
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Figure 5. - Load distribution due to pitching oscillations about y-axis 
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Figure 5. - Concluded. Load distribution due to pitching oscillations about y-aocis. 
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